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We consider a generalization of the Dicke model. It de- 
scribes localized, physically separated, saturable excitations, 
such as excitons bound on impurities, coupled to a single long- 
lived mode of an optical cavity. We consider the thermal 
equilibrium of the model at a fixed total number of excitons 
and photons. We find a phase in which both the cavity field 
and the excitonic polarization are coherent. This phase corre- 
sponds to a Bose condensate of cavity polaritons, generalized 
to allow for the fermionic internal structure of the excitons. 
It is separated from the normal state by an unusual reentrant 
phase boundary. We calculate the excitation energies of the 
model, and hence the optical absorption spectra of the cavity. 
In the condensed phase the absorption spectrum is gapped. 
The presence of a gap distinguishes the polariton condensate 
from the normal state and from a conventional laser, even 
when the inhomogeneous linewidth of the excitons is so large 
that there is no observable polariton splitting in the normal 
state. 
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I. INTRODUCTION 

In the strong-coupling regime for matter and light, ra- 
diative decay of a material excitation gives way to cou- 
pled oscillations of the polarization of the matter and of 
the electromagnetic field. The quasiparticles correspond- 
ing to such coupled modes are known as polaritonsEJ. 
The classic realization of polaritons is excitons in a bulk 
semiconductor coupled to photons inJree space, as dis- 
cussed many years ago by HopfieldEl. In this exam- 
ple, wavevector conservation ensures that each exciton 
is coupled only to a single mode of the electromagnetic 
field, leading to the formation of polaritons which are su- 
perpositions of a single exciton and photon. Recently, 
there has been a lot of interest in polaritons formed 
from photons confined in cavities: such cavity polari- 
tons have now, been observed for confined photons cou- 
pled i.o atomsu, to two-dimensional excitons in quantum 
wellsQ, to bulk .egcitonsEl, to excitons in films of organic 
semiconductorsBQ, and to charged exciton complexes!! 

Since polaritons are photons coupled to other excita- 
tions, they axe bosons, and so are pandidates for Bose 
condensation!! Recent observationsEErtJ of bosonic be- 
havior for cavity polaritons have renewed interest in this 
idea. 

However, there is a conceptual difficulty with a Bose 
condensate of polaritons: while polaritons are usually 
considered in the low-excitation linear regime, Bose con- 



densates are stabilized by nonlinearitiesEl. For cavity 
polaritons, there is also the following more practical dif- 
ficulty. Bose condensates are characterized by coherence, 
and in a polariton condensate this coherence will appear 
in the photons. Given this, how is a polariton condensate 
distinct, conceptually and observationally, from a laser? 

In this paper we address these problems by developing 
a theory of polariton condensation in the Dicke modclEl 
This nonlinear model of confined photons coupled to mat- 
ter is one of the basic models of laser physics. It allows 
us to go beyond the conventional linear-response concept 
of a polariton, including effects due to finite excitations 
of the matter in the cavity. In the language of semicon- 
ductors, it includes a "saturation" or "band-filling" non- 
linearity, produced by the fermionic internal structure of 
the excitons. 

Polaritons are not conserved particles, so there is ul- 
timately no equilibrium condensate. We may, however, 
treat polaritons as conserved particles if their lifetime is 
much longer than the time required to achieve thermal 
equilibrium at a fixed polariton number. We will study 
this quasi-equilibrium regime, since it. is in this regime 
that Bose condensation is well-definecfi. 

In section ||, we introduce the model, and explain how 
the concept of a polariton can be generalized to allow 
for the nonlinearity of the model. We then present, in 



sections [II 



and IV, a simple variational technique for 
calculating the ground state of the model at a fixed den- 
sity of polaritons. In section ^ we investigate the ther- 
modynamics of the model using an alternative technique 
based on functional integrals. This technique demon- 
strates that the variational approach is essentially exact, 
and allows us to consider finite temperatures. In sec- 
tion VI , we use the expressions derived by the functional 



integral method to stud y th e phase diagram for conden- 
sation, while in section VII, we use these expressions to 
calculate the excitation spectra of the model. These ex- 
citation spectra provide a physical picture of the transi- 
tion to the condensed state, and determine the absorp- 
tion spectrum of the cavity. Finally, in section VIII we 
discuss our conclusions. 

The functional integral approach to the thermodynam- 
ics of p«r model has already been the subject of a brief 
reportliZl. We extend that earlier report to allow for a dis- 
tribution of the energies of the electronic excitations, i.e. 
inhomogeneous broadening, which is significant in many 
potential realizations of the polariton condensate. 
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II. MODEL 



The Dicke modeOll3 consists of a set of N two-level os- 
cillators coupled to a single mode of the electromagnetic 
field by the dipole interaction. The two-level oscillators 
do not interact with one another, except through their 
common coupling to the electromagnetic field. We gener- 
alize the original Dicke model to include an energy distri- 
bution of the two- level oscillators. Making the rotating- 
wave approximation(see e.g. Refs. |i^ , |l9| ), we consider 
the Hamiltonian 

H = J2^^ (tfb-<Ja)+uj c tf4> + H', (1) 

H' = -^^(bW + vVb) . 

Here the two-level oscillators are indexed by the variable 
n, which is summed over. We use a fermionic represen- 
tation for the two- level oscillators, describing each one 
in terms of a pair of fermions with annihilation opera- 
tors a and b. For brevity we suppress the index n on 
the fermionic operators. The fermions are subject to the 
single-occupancy constraint 



tfb + a<a = 1 



(2) 



on each site, ip is the annihilation operator for the cavity 
mode, E g (n) is the energy of the n th two-level oscillator, 
and g is the strength of the dipole coupling. 

The Hamiltonian (Q) is a simpk^ppdel of a three- 
dimensional cavity (photonic dot)Ej£!l, containing lo- 
calized, physically separated electronic excitations. Al- 
though simplified, it is a useful starting point for many 
systems. For example, each of the two-level oscillators 
could describe the presence or absence of a localized ex- 
citon in a given eigenstate of the disorder potential in a 
disordered quantum well, on a given molecule in an or- 
ganic film, or trapped on a particular impurity. The re- 
striction to singly occupied states models the hard-core 
repulsion produced by the fermionic structure of such 
excitations. It describes spinless excitations which are 
localized on the scale of the (exciton) Bohr radius. It 
is straightforward to generalize our calculations to allow 
for a finite number of excitations on each site, i.e. traps 
which are bigger than the Bohr radius, and so can hold 
several excitons. 

H' couples the photons to excitations of the two-level 
oscillators, created by the operator S+ = -j=Yl^ a. If 

Eg {n) = Eg then the excited states which are created by 
S + from the vacuum are eigenstates of the bare Hamil- 
tonian H. If furthermore N is large and the two-level 
oscillators are near to their ground state, 



i^(6t6- a t a )«-iV/2, 



(3) 



cillatorsQ. Polaritons are usually presented as the eigen- 
states of such a model. 

Away from the low-excitation limit (||), S+ is not a 
bosonic creation operator, and the conventional descrip- 
tion of polaritons breaks down. To go beyond the low- 
excitation limit, we generalize the concept of a polariton 
to be the quantum of excitation of the coupled matter- 
light system. The polariton number is then the total 
number of photons and excited two-level oscillators, 



N pol = L + N/2 



(4) 



a 1 a + N/2, 



which is a conserved quantity for the model ([j]). de- 
fines the operator L, which we refer to as the excitation 
number. We define a corresponding excitation density 
p ex — (L)/N, which is the total number of photons and 
electronic excitations, per two-level oscillator, minus one- 
half. Since the numbers of photons and electronic excita- 
tions are positive, the lowest excitation density is —0.5. 
Since the number of electronic excitations is always less 
than N, the electronic contribution to p ex is always less 
than 0.5. 

The thermal equilibrium of the Dicke model, in the 
absence of an externally created population of polari- 
tons, has been studied extensively since the pioneering 
exact solution of Hepp and LiebE3. These authors showed 
that, even in the absence of external excitation, the Dicke 
model has a phase transition to a Bose condensed state. 
Such an equilibrium condensat&is a static, coherent state 
of photons: it is a ferroelectricca. Here we are interested 
in the thermal equilibrium of a population of polaritons: 
the quasi-equilibrium problem posed by (^) at a fixed ex- 
citation L. The quasi-equilibrium condensate which we 
find in this regime is a time- varying generalization of the 
ferroelectric state discovered by Hepp and Lieb. 



III. VARIATIONAL APPROACH 

We can write down a variational state which describes 
the polariton condensate by noting that Bose conden- 
sates are described by coherent states. This produces 
a variational wavefunction closely related to the BCS 
wavefunction used to describe superconductors and ex- 
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ensates. As has been stressed by Comte and 
, this class of wavefunction can describe an ex- 
citon condensate in both the low and high density limits. 
It thus permits a smooth interpolation from low densities, 
where the excitons are simple bosons, to high densities, 
where their fermionic internal structure is revealed. In a 



then S+ is approximately a bosonic creation operator, 
and the Hamiltonian ([!]) becomes two coupled boson os- 



1 This mapping may be formalizedta using the Holstein- 
Primakoff transformation. 
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similar way, it allows us to explore the polariton conden- 
sate beyond the low-density regime in which polaritons 
are usually considered. 

The constituents of our proposed condensate are the 
excitations of the whole system, counted by L. In gen- 
eral, such an excitation is a superposition of an excitation 
of the cavity mode and an excitation of the electronic 
states. Thus we take for our trial wave-function a coher- 
ent state of such a superposition: 



| A, id) = e 



\vac) 



(5) 



where the state \vac) has a single fcrmion in the lower 
state of each two-level oscillator. The state (||) has a 
finite polarization of the electronic excitations as well as 
a finite amplitude for the cavity field. A and w n are the 
variational parameters. Expanding the exponential, (j^) 
explicitly becomes a superposition of a coherent state of 
photons and a BCS state of the fermions, 



|A, u, v) 



'[]Mt+ Un e*flt)|0). 



(6) 



Here \,u n ,v n , and <j) n are the variational parameters, 
and |0) denotes the vacuum state with no fermions in 
any of the levels. By construction, this variational state 
obeys the single-occupancy constraints (||). We fix the 
overall phase of the condensate by choosing A to be real. 
The 4> n have been explicitly introduced to make the u 
and v real. They are the phase differences between the 
cavity field and the polarizations of the electronic states. 

To find the ground state of (Q) at fixed excitation num- 
ber we minimize 



(H - HexL) = ti c A 2 + e n (v 2 n - u 2 n ) 



+2^=\u n v n cos(0 n ), 
v N 



(7) 



E g (n) - fi ex 



with respect to the variational parameters, subject to the 
normalization conditions u 2 + v 2 = 1. 

Although the overall phase of the condensate is arbi- 
trary, the relative phases 4> n are not: there is only one 
order parameter. The relative phases <p n are fixed by 
the last term in ([?]), the dipole coupling. This term en- 
sures that all the two-level oscillators which have a finite 
dipole moment (u„ ^0,1) are mutually coherent, 4> n = 4>, 
when the energy is minimized. It is the dipole interac- 
tion which is responsible_£or Bose condensation, and its 
accompanying coherencetJ, in the present system. 

Setting 4> n = and defining an intensive A by rescaling 
A — » X\^N, the condensate parameters are given by the 
real solutions with Xu n v n < to 



w c A + -y u n v n = 0, 

n 

2e n u„v n - g\{v 2 n - u 2 ) = 0. 



(8) 



\i ex was introduced as a Lagrange multiplier constrain- 
ing the excitation number. It is the chemical potential 
for our coupled modes, and is related implicitly to the 
excitation density by 



= a 2 + — vi -ui. 



a) a 



(9) 



Eliminating u n and v n from (||) and (^) we can rewrite 
these expressions as 



^A V J_ 



Pe X ~ A 2 - — J2 



2N^\E n \' 
where we define 

E n = sign(e n )v / i 2 Ti 2 |A| 2 - 



(10) 

(11) 

(12) 



( |l0| ) is analogous to the BCS gap equation, with an order 
parameter A. 



IV. ZERO TEMPERATURE PROPERTIES 



To investigate the expressions (|10|-|12[), we replace the 
summations over sites with an integral over the energy 
distribution of the two-level oscillators. We take this dis- 
tribution to be a Gaussian with mean Eq and variance 
ag. The remaining parameters in our quasi-equilibrium 
problem are then the excitation density, p ex , and the 
dimensionless detuning between the energy of the cavity 
mode and the center of the exciton line, A = (uj c — Eq) / g. 

For a Gaussian density of states, the summation on the 
right of (|Io| ) diverges as A — > 0, and approaches zero as 
A — > oo. Thus for any /i ex < uj c there is always a con- 
densed solution, A 7^ 0, to (10): the system is condensed 
at arbitrarily small excitation densities. This behavior is 
produced by the tails of the Gaussian distribution. Be- 
cause of these tails, we have excitons at arbitrarily low 
energies, and hence also bound exciton-photon states at 
arbitrarily low energies. It is impossible to populate just 
the excitons, because no matter how small [i ex is, there 
is always a bound state involving photons below it. Wc 
expect that if the density of states has a lower cut-off, 
and is continuous at this cut-off, there would be a finite 
critical [i ex below which there is no condensed solution 
to©. 

Let us investigate the dependence of fi ex on p ex in the 
absence of inhomogeneous broadening, a — 0. At low 
densities, p ex —0.5, p ex can be obtained from (|lC|). 
Expanding this expression for small A and comparing the 
leading terms, we find that fi ex is given by the conven- 



tional linear-response polariton energy, \i e 
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FIG. 1. Right panel: dependence of the chemical potential 
on excitation density for detunings A = 0, 1 and 3 and vari- 
ances a = 0,0.5 and 1. Left panel: absorption spectrum for 
a microcavity at p ex = — 0.5 and T — for a — 0.5 and the 
same three detunings. 



i [(uj c + Eg) — gyj A 2 + 4] . At finite densities we calculate 
/j, ex numerically, by solving ( |To| ) and ( p"T| ) to determine 
Pex(f^ex)- The results are plotted in the right panel of 
Fig. |l|, for A = 0, 1 and 3. At low densities we are de- 
scribing a condensate of conventional polaritons, and so 
have n ex = i^LPB ■ As the density is increased the exciton 
states saturate, forcing the excitations to become more 
photon-like. Thus the chemical potential approaches co c 
at high densities. For A > 2 the separation between 
the exciton-like and photon-like excitations persists to 
Pex — 0.5, where the exciton states are completely satu- 
rated. This results in a discontinuity in p ex at this point, 
since no further excitation can be added to the exciton 
states. 

The dependence of jj. ex on p ex in the inhomogeneously 
broadened case is also illustrated in the right panel of 
Fig. |l|. It is qualitatively rather similar to the homoge- 
neous case. Instead of the finite intercept of the homo- 
geneous case we now have jji ex ^ -oo as p ex — ► —0.5. 
This behavior is again caused by the tails of the Gaus- 
sian distribution. To demonstrate how p, ex approaches 
the conventional polariton energy -Elpb in the homoge- 
neous, low-density limit, we compare the behavior of p ex 
with the density of states for the linear-response excita- 
tions of the empty(p ea; = —0.5) cavity. This density of 
states is the optical absorption spectrum of the cavity, 
and is plotted in the left panel of Fig. [I] for a = 0.5 and 
A = 0, 1 an d 3. We will describe how it is calculated in 
section VII . At very low densities, \x ex lies in the tails of 
the exciton distribution. With increasing density, these 
states quickly saturate, producing a sharp rise in p, ex . As 
fi ex reaches the polariton peak, the sharp rise in the den- 
sity of states for the coupled modes produces a kink in 
the chemical potential. In the homogeneous limit, this 
kink moves to zero density and corresponds to the usual 
polariton energy. Since the density of states at this point 
is infinite in the homogeneous limit, these polaritons are 




(E-E )/g 

FIG. 2. Occupation of the two-level oscillators at zero tem- 
perature as a function of energy E for A = 3, T = 0, a = 0.5 
and densities p ex = —0.4, —0.2, 0, 0.2, 0.4, 0.6(dot-dashed 
curves, increasing from left to right) and p ex = 100(dotted 
curve) . The shaded region shows the Gaussian distribution of 
oscillator energies used. 



simple bosons. 

Figure || shows the occupation of the two-level oscil- 
lators in the polariton condensate, for A = 3, a = 0.5 
and various densities. The occupation number of the n th 
two-level oscillator is 



As is clear from the figure, this is a Fermi step broadened 
by the interaction with the photons, just as the elec- 
tronic distribution in a BCS superconductor is a Fermi 
step broadened by the pairing interaction. The states 
in the broadened region of the step have a finite dipole 
moment and are involved in the condensate. The Fermi 
step moves up through the exciton line as the excita- 
tion is increased from p ex — — 0.5 and the low- lying elec- 
tronic states saturate. At very large densities there are 
a large number of photons, and the Fermi step is almost 
completely flat: rather than the electronic system com- 
pletely saturating in the high density limit, it approaches 
half filling. This is because the half-filled state maximizes 
the polarization of the electronic states and hence min- 
imizes the dipole interaction between the excitons and 
the macroscopically occupied cavity mode. 

Careful inspection of Fig. || reveals that the broadening 
of the Fermi step produced by the photons does not in- 
crease monotonically with density. This corresponds to 
a non-monotonic dependence of the field amplitude, A, 
on density. This dependence is illustrated in Fig. |[ The 
field amplitude is related to the electronic polarization 
by the first of the equations (||). It is proportional to the 
electronic polarization and inversely proportional to the 
separation between the chemical potential and the cavity 
mode. The electronic polarization depends on the den- 
sity of states in the vicinity of the chemical potential(Fig. 
||); the peak in the density of states at the center of the 
exciton line produces the peak in Fig. 0. 
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FIG. 3. The order parameter A as a function of density, for 
a — 0.5 and A = 0, 1 and 3. A 2 is the photon number per 
two-level oscillator in the condensed state. 



V. LARGE-N EXPANSION 



The variational approach of sections II] and IV be- 
comes exact in the thermodynamic limit N — > oo. Physi- 
cally, this is because it corresponds to a mean-field treat- 
ment of the interaction between electronic excitations. 
This interaction, between a large number (iV) of elec- 
tronic excitations, is mediated by a small number (one) 
of cavity modes. In a mean-field treatment of this inter- 
action, each electronic excitation is coupled to the aver- 
age field produced in the cavity by the other electronic 
excitations. This becomes exact when there are a large 
number of electronic excitations contributing to a small 
number of field modes, since the fluctuations of the field 
are then negligible. 

In this section, we develop a mean-field theory for the 
thermodynamics of the model ([!]) from the functional- 
integral representation of the partition function. In this 
representation, the partition function can be rigorously 
evaluated, for large TV, using a saddle-point analysis^. 
From such an analysis, we derive finite-temperature gen- 
eralizations of the variational expressions (fL0|-|l2|), thus 
demonstrating that they are rigorous in the limit of large 
N. 

The functional . in t eg ral techniques used here have pre- 
viously been usedc3'E3 to calculate the partition function 
and excitation energies in the absence of a constraint 
on the polariton number of a simplification of the Dickc 
model. While the Hamiltonian of the model discussed in 
Refs. |6|,|7| is given by (|l|) with E g (n) = E g , the local 
constraints prohibiting two fermions on the same site, 
(||), are replaced with a global constraint. In contrast, 
we retain (0) as local constraints, as well as including 
a distribution of E g and a constraint on the polariton 
number. 

As in sections [II and |y|, we work in a grand-canonical 
ensemble, using a chemical potential /i ex to constrain the 
excitation number. We consider the partition function 



associated with this ensemble 

Q = Tre-^ H -^ L \ 

The coherent-state functional-integral formalism allows 
us to express Q, for the model ([!]), as the constrained 
functional integral 



Q = J V^]J[DT] n S(fj n r, n - l)]t 



with the action 

S = drip{d T + u c )ip + ^ f\ n M n r] n . 

We have introduced a Nambu spinor 

( bn 

for each two-level oscillator. The matrix M n is 



M _ ( d T +i n gtp/VN 
\9$/VN d r -e n 



Rescaling the boson field ij> — > vNip and transferring 
the fermionic integrals into the action gives 



Q = J V^\J\e 



NS 



with an effective action 

rP 



f 1 

Sett = / dTlp(d T + U! c )lp - — ^ Sf, n , 
n 

S{, n = In [ Vrj n S(fj n ij n - l)e ^° nnPnV " , 



(13) 



in which the P n are the matrix operators M n after rescal- 
ing the boson field, and J denotes the trivial Jacobian 
arising from this rescaling. 



A. Mean-field equation 

For large N, the dominant contribution to the parti- 
tion function Q comes from those functions i/jo(t) which 
minimize the action S c ff. Such functions obey the Euler- 
Lagrange equation. For the action (|l3|), this takes the 
form 



(d T +lj c )?Po(t) = — 



i/!(t)=V>o(t) 



N t—' dtp 



-jj^2( a n{r)b n (T)) 7 



(14) 



where the right-hand side of this expression is the polar- 
ization of the two-level oscillators in thermal equilibrium 
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driven by an external field %Pq(t). This polarization ap=. 
pears because the field ^o(t) modifies the eigenstatesEj 
of the electronic system. A thermal population of these 
new eigenstates can correspond to a finite polarization of 
the original fermions. Equation ( fl4| ) is a self-consistency 
condition: the cavity field is driven by the polarization of 
the fermions, which itself arises from the renormalization 
of the fermions produced by the photons. 

Assuming that the self-consistent field ipo( T ) 1S inde- 
pendent of t, we can calculate the polarization term on 
the right of (111) by making a Bogolubov transformation 



cos(0)e^ - sin(0) 
sin(0) cos(0)e"^ 



tin 

In 



(15) 



from the b n and a n fermions to new fermions 5 n and "f n . 
This transformation diagonalizes P n when 4> = arg A and 
tan 20 = g\\\/e n . The S n and j n quasiparticles then have 
energies ±E n respectively, with E n defined by equation 
(|l2|). Since ( |l5| ) is a rotation in 77 space, it preserves the 
single occupancy constraints. Thermally populating the 
new fermions in accordance with the single occupancy 
constraint we have 

(a n b n ) = ie 40 sin(20)(J„5„ - %y n ) 



-e^ sin(20) tsnh(J3E n ), 



and (114) becomes 



A =§Ei-^h(^). 



(16) 



Equation ( |l6| ) is the finite-temperature generalization 
of the variational result (|To|). This generalization is 
rather straightforward: we have just acquired tanh(/3i?) 
factors describing the thermal occupation of the two-level 
oscillators. 

If we remove the constraint on the polariton number, 
by setting [i ex = 0, and set E g (n) = E g , then ( |l6| ) is 
the form originally derived by Hepp and LiebE3 for the 
unconstrained equilibrium of the Dicke model. In that 
problem, the existence of a condensate requires 



COr.Ea 



< 1, 



(17) 



since otherwise (|l6|), with \x ex — and E g (n) = E g , has 
only the trivial solution A = 0. However, it is shown 
in Refs. ^ and |3(] that the A 2 terms of the minimal- 
coupling Hamiltonian, neglected in the model (^) , modify 
the inequality (|l7]) in a way which is inconsistent with the 
Thomas-Kuhn-Reich sum rule. This sum rule requires 
KEg/g 2 > 1, where k is the coupling constant for the A 2 
term, while the modified inequality (|l7|) reads 



2k)E 



< 1. 



r 



(18) 



Since this inequality cannot be satisfied, the phase tran- 
sition in the unconstrained case is an unphysical artifact 
of the model ([j]). However, we do not believe that the 
A 2 terms prevent condensation in the constrained case, 
because the inequality corresponding to ( |T§| ) will be 
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< 1, 



and the parameter \i ex is not restricted by the sum rule. 



B. Effect of fluctuations 

Let us now consider the effect of small fluctuations 
5ip(r) around the mean-field solution. Expanding S c s 
to second order in a functional Taylor series around the 
mean-field solution we have 



Q 



-NS 



T>(Si/j) I J\e~ NS2 1 5 ^' 5 ^'] . 



(19) 



Here 5*o is the action evaluated on the extremal trajectory 
and S*2 is the quadratic action from the second order term 
in the Taylor series. S2 is the effective action for small 
fluctuations of the electromagnetic field. The kernel of 
S2, G^ 1 , is the inverse of the thermal Green's function 
for the photons. 

The integral over fluctuations in ([n]) contributes a 
term 



1 

N 



In det Q 



to the free energy density. Since the mean-field solution 
should be a minimum of the action, the eigenvalues of 
Q^ 1 should be positive. Then In det C/ -1 is finite as N — > 
co, there is no fluctuation contribution to the free energy 
density in this limit, and the mean-field theory becomes 
exact. 



C. Effective action for fluctuations 

However, we have yet to check whether the solutions to 
( |l6| ) are actually minima of the action or merely extrema, 
i.e. whether the mean-field solutions are stable against 
fluctuations. To check this, we will need the effective 
action £2, which we derive in this section. 

To obtain S2, we calculate the (functional) second 
derivatives of S c g, and evaluate them on the extrema 
4>{t) = iPo{t) = A. In the frequency representation, the 
components of S2 are 



d 2 S cS 



(20) 



9 -y. 

N ^ 



:)«)) 



dr, 



G 



and 



d 2 S cff 



- = -f3g 2 5(to' +co) 



(21) 



f eiUT «^(r)^(0)) 

n ^° 



lj and a/ denote bosonic Matsubara frequencies, cr+ = 
6^a„ is the polarization operator for the n th two-level 
oscillator, and the integrands are the susceptibilities of 
the two-level oscillators in the self-consistent field A. 

( po|) and ( plf ) describe coupled fluctuations of the cav- 
ity field and the electronic polarization. They .axe anal- 
ogous to the Dyson-Gor'kov-Beliaev equationsBj of the 
theory of superconductors and weakly-interacting Bose 
gases. Because we are working with a condensed system 
there is an anomalous contribution to 52, (pl|), from fluc- 
tuations which do not conserve the number of excitations 
above the condensate, i.e. those in which an excitation 
enters or leaves the condensate. 



To calculate the susceptibilities which appear in (20), 
we rewrite them in terms of the renormalized two- 
level oscillators using the transformation (fl5|), trans- 
form to the Schrodinger representation, and take thermal 
and quantum-mechanical averages over the renormalized 
eigenstates. This gives 



S2[Sip, Sip] 

g- 1 



0£ W) W-»))Q- X ( 



K x K 2 
K* K\ 



(22) 



t 
N 

ie n oj - 2el 



Ki = iu + o> e + — 



— tanh (J3E n ) 

tin 



5 2 |A| 5 



4\2 



g*X 
N 



E 



E n {uj^+AEl 



$u>u n \\\ 2 g 2 
- tanh {j3E n 



(3 



sech" (/?£„) . 



Note that, for the condensed state, Q~ x takes different 
forms at to — and at finite lj. This is because thermal 
fluctuations at co = include both fluctuations of the 
order parameter and quasiparticle excitationsEil. Only 
the latter appear at finite lo. In the normal state, A = 0, 
the effective action simplifies to 



5 2 = ^E^H 



ILU + LJ C 



1 v 

N ^ 



1 • 

9i 



2ie. 



— tanh (f3e n ) 



(23) 



D. Nature of the extrema 



We now use the expressions (p!2|-|23[) to investigate the 
nature of the extrema when E g (n) — E g . 

Considering first a condensed solution, A ^ 0, we use 
the extremal equation ( |l6| ) to eliminate -g- tanh ((3E n ) 

from the matrix Q~ x . The eigenvalues of the resulting 
matrix are all strictly positive provided that lj c > 0, ex- 
cept for a single zero eigenvalue at uj n = 0. From jl^ ) 
we see that the condensed solutions always have tD c > 0. 
Thus we conclude that, at a condensed solution, the ac- 
tion has a minimum in all but one direction, and is locally 
flat in this one direction. 

We show in Appendix |A| that the single zero eigen- 
value describes a change in the overall phase of the con- 
densate. It is the Goldstone mode corresponding to the 
broken gauge symmetry of the condensate. Because we 
are considering a broken symmetry state, we should not 
integrate over these fluctuations when calculating the 
partition function. Since the other eigenvalues of Q~ l 
are always positive for the condensed solutions, these so- 
lutions are stable against physical fluctuations, and the 
mean-field theory is exact|^|. 

Turning now to the normal solution, A = 0, we find 
from (123) that this is a minimum of the action unless 



£> c < — tanh(/3e) 



(24) 



This is just the condition for the extremal equation jl^ ) 
to have a condensed solution. Thus we have the usual sce- 
nario of a continuous phase transition: there is a phase 
boundary (24), at which the normal state becomes un- 
stable and a stable, condensed solution appears. 



E. Density equation 

As well as the mean- field equation ([Trf), we need the 
equation relating the density p ex to the corresponding 
chemical potential fi ex . This is obtained from the parti- 
tion function in the standard way, 



d 



(3N dfl e 



InQ. 



(25) 



2 In Appendix [a|, we give a formal demonstration that the 
zero mode does not contribute to the free energy density as 
N — » oo, so that the presence of th e zer o mode does not 
invalidate the discussion of subsection VB. 
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FIG. 4. Phase boundaries for A = 0(left panel), A = l(center panel) and A = 3(right panel), and variances a = 0(solid 
lines), a — 0.5(dotted lines) and a = l(dashed lines). For A = 3, a = 1 the upper branch of the phase boundary lies off the 
scale, while for A = 3, a = 0.5 the lower branch is indistinguishable from the homogeneous case. 



The asymptotic form for the partition function is Q ~ 
e~ NS °, where So is the minimal action. Inserting this 
asymptotic form in ( p5| ) gives, for the solution ipa( T ) = 



Pe 



|A| S 



2N ^ EL 



tanh (@E n ) , 



(26) 



which is the generalization of ( pi] ) to finite temperatures. 

The first term in (26) is the contribution to the excita- 
tion density from the macroscopic electromagnetic field, 
while the second term is the contribution from the ther- 
mal population of renormalized electronic excitations. In 
the absence of a macroscopic electromagnetic field, A = 0, 
both the photon contribution and the renormalization of 
the electronic excitations disappear. The expression (Pq) 
is then the familiar form for the excitation of a set of 
two-level oscillators. 



VI. PHASE DIAGRAM 

From (p4| ) and (^6|) we have the critical temperature 
for condensation, as a function of the excitation density, 
in the homogeneous model: 



Peg = 



4 tanh 1 (2p c 



A ± VA 2 - 8p e 



(27) 



Note that the transition temperature depends logarith- 
mically on the density, and its scale is set by the inter- 
action strength g. This is in contrast with a model of 
propagating, weakly-interacting bosons, where the tran- 
sition temperature varies as a power law of the density 
and its scale is set by the mass of the bosons. 

At low densities, (2?j) is the phase boundary separat- 
ing a population of electronic excitations with energy Eq 
from a population of conventional polaritons with en- 
ergy -Elpb- To see this, note that such a transition 



would occur when the chemical potential for the elec- 
tronic excitations reaches -Elpb , corresponding to a den- 
sity p ex + 0.5 « e -/3o(£o--BLPB) ; w hich is the low-density 
limit of (§7|). 

For the inhomogeneous model, we calculate the phase 
boundary numerically, assuming the same Gaussian dis- 
tribution of energies as in section IV. We obtain the 



critical chemical potential for condensation, /i c (/3 c ), by 
demandingthat ([!(]) have a repeated root A = 0, and 
then use (Era) to obtain the critical density p c (f3 c ). 

In Fig. ^pwe plot the homogeneous phase boundaries 
( p7| ) , along with numerical results for the inhomogeneous 
model with a = 0.5 and 1. On resonance, A = 0, the 
transition temperature increases monotonically with den- 
sity. The system is always condensed for p ex > 0, because 
to exceed this density would require a chemical potential 
above the center of the energy distribution of the elec- 
tronic excitations, and hence above the bosonic cavity 
mode. While for A < 0(not illustrated) the phase bound- 
ary is qualitatively unchanged from the resonant case, for 
A > we find reentrant behavior. This behavior is the 
result of the saturable nature of the electronic states. It 
can be understood by considering the limits p ex — > ±0.5 
when A 3> 0. Near the p ex — —0.5 limit, the normal 
state consists of a small number of electronic excitations, 
weakly interacting with each other through the cavity 
mode. They condense when their density exceeds a crit- 
ical value set by the strength of their interaction, which 
is determined by A and g. Near the p ex — 0.5 limit, the 
electronic system is constrained to be fully occupied, and 
the normal state consists of a small number of holes in 
an otherwise completely excited electronic system. These 
holes again interact through the cavity mode, and so the 
transition occurs when the density of holes, 0.5 — p ex , ex- 
ceeds a critical value. For A — > oo the critical densities 
of holes and excitons are identical, so the phase diagram 
is symmetric about p = 0. For finite A, the interaction 
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a a 
FIG. 5. Dependence of the critical density on the inhomo- 
geneous broadening er, for A = 0(left panel) and A = 1 (right 
panel), and /3 = l(top curve), 3, 5, 7, 9, 11, 13(lowest curve). 



is stronger for the holes than for excitons, since they are 
nearer in energy to the cavity mode, and so the phase 
boundary becomes skewed to the forms shown. 

At temperatures which are high compared with the in- 
homogeneous broadening erg, thermal fluctuations domi- 
nate over the inhomogeneous broadening. Thus at these 
temperatures the inhomogeneous broadening has little 
effect, as can be seen in Fig. [|. However, at low tem- 
peratures the inhomogeneous broadening suppresses con- 
densation by increasing the energy separation between 
the electronic excitations and the photons, collapsing the 
phase boundaries towards p ex = 0. The effects of inho- 
mogeneous broadening are further illustrated in Fig. |^, 
which shows the dependence of the critical density on a 
at various temperatures for detunings A = and A = 1. 



VII. EXCITATION ENERGIES 



In this section, we use (|22|) and (|23|) to study the ex- 
citation spectrum of the quasi-equilibrium states of the 
model (Q). The excitation spectra we calculate explain 
the form of the phase diagrams in Fig. [|. The excitation 
spectra of the two quasi-equilibrium states are different 
from each other, and also from the excitation spectrum 
of a conventional laser. Since the excitation spectrum 
is directly related to the optical absorption spectrum of 
the cavity, which is an experimentally accessible quan- 
tity, these spectra offer a clear experimental signature of 
polariton condensation. 

The matrix Q" 1 , given by (p2|), is the inverse of the 
thermal Green's function for the photons. We use the 
standard relationsOO between thermal Green's func- 
tions, retarded Green's functions, and excitation spectra, 
to extract the latter from (fi2|). 



A. Homogeneous model 

We begin with the normal state of the homogeneous 
model. The inverse of the normal state Green's function 



contained in (|23|) can be written as a sum of simple poles 



iuj n + E+ iu> n + E- 



(28) 



The structure of this Green's function is clear: we have 
two excitations, with quasiparticle energies 



E± + Mcx = [K + E g ) ± g^A 2 - 8p ex ]/2, 

and corresponding weights 

C± =±(2i-E±)/(E_-E + ). 

These normal-state excitations are polaritons in the gen- 
eral sense of Hopfieldo: coupled modes involving the lin- 
ear response of the electronic system around its equilib- 
rium state. The gap in the spectrum is increased over the 
bare detuning A owing to the dipole coupling between the 
excitons and the cavity mode. The presence of excitation 
in the ground state, either driven by finite temperatures 
or by finite /i C x, causes the two polariton branches to 
attract. This attraction is due to the decrease in the po- 
larizability of the electronic states as their population in- 
creases and saturation occurs. It can also be understood 
in terms of an angular momentum representationE3 for 
the collective states of the electronic system. In such a 
representation, the excitation of the electronic states cor- 
responds to the z component of an angular momentum, 
while their polarization corresponds to the raising oper- 
ator S+. Thus the polarizability of the electronic states 
is maximized at (S z ) — —N/2. 

Since condensation is a phase transition, we expect 
a qualitatively different excitation spectrum in the con- 
densed state. From ( p2| ) and ([!(]) , we find for the leading 
component of the matrix thermal Green's function 



£j c {uj 2 + 2g 2 \\\ 2 )- iu(uj 2 + AE 2 + 2Q c e) 



(iuj n ) 2 (iu; n + £)(iuj n -£)(! + 5 Un a) 



(29) 



with £ = y/{£j c + 2e) 2 + Ag 2 \X\ 2 . The interpretation of 
this Green's function is complicated because, as we have 
already mentioned, it describes both quasiparticle exci- 
tations and fluctuations of the order parameter. To rig- 
orously obtain the quasiparticle spectrum, we should ex- 
tract the contribution to (^9|) from the order parameter 
fluctuations, and analytically continue the remainder to 
obtain the retarded Green's function. Rather than fol- 
low such a procedure, we propose the following physically 
appealing if mathematically naive interpretation of (|2^) : 
the Kronecker delta and the (iuj n ) 2 terms in the denom- 
inator describe the condensate response, leaving quasi- 
particle excitations at energies ±£. The (iu> n ) 2 is clearly 
associated with the phase mode of the condensate, dis- 
cussed in Appendix while the Kronecker delta is re- 
lated to number fluctuations of the condensates. The ex- 
citations at energy ±£ are coupled exciton-photon modes 
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Pex Pex 

FIG. 6. Excitation energies and chemical potentials as a 
function of density for the homogeneous model at A = 
(left panel) and A = 2 (right panel), both with gf3 = 2. Thin 
solid lines: normal state excitation energies. Thick solid lines: 
condensed state excitation energies. Dashed lines: normal 
state chemical potential. Dot-dashed lines: condensed state 
chemical potential. The shading marks the condensed region 
for this (5. 



of Fig. |6| the chemical potential crosses the lower polari- 
ton branch at p cx = without the condensate appearing. 
It is not until the chemical potential crosses the upper 
polariton branch that the transition occurs. This can 
be understood by considering the signs of the quasipar- 
ticle weights C±. A positive quasiparticle weight corre- 
sponds to absorption of an external field (a particle- like 
transition), whereas a negative quasiparticle weight cor- 
responds to gain(a hole-like transition). For p c ^ > 0, 
the lower polariton branch has a negative weight: it has 
become hole-like, and must be below the chemical poten- 
tial for stability. At the transition it is now this lower 
branch which joins to the "pair forming" excitation of 
the condensate, while the upper branch joins to the phase 
mode and the "pair breaking" excitation appears above 
the phase mode. 

B. Inhomogeneous model 



in the presence of the macroscopic electromagnetic field 
of the condensate. £ is analogous to the pair breaking 
energy in a superconductor: it is the energy required to 
extract an exciton-photon complex from the condensate. 
Note that if we remove the photon contribution to this 
energy, byj-setting uj c = 0, then £ becomes the familiar 
expressioncS for the energy of an electron-hole pair in the 
presence of a classical electromagnetic field at frequency 

In Fig. ||, we illustrate the evolution of the excitation 
energies of the microcavity with increasing density. To 
explain the relationship between the excitation energies 
and the phase diagram, we also plot the chemical poten- 
tials for the normal and condensed states on this figure. 
The left panel of this figure should be compared to the 
g/3 = 2 line of the corresponding phase diagram, which 
is the left panel of Fig. |J. When A = and p cyi = — 0.5 
the system is in the normal state. Increasing p ex pop- 
ulates the electronic excitations, increasing the chemical 
potential and decreasing the polariton splitting. Even- 
tually the chemical potential crosses the lower polariton 
branch from below and the system condenses. At the 
critical density, the lower polariton branch joins to the 
phase mode at the chemical potential, the upper branch 
joins to the "pair breaking" excitation, and an excitation 
appears below the chemical potential. This latter exci- 
tation has zero weight at the transition. It corresponds 
to an excited state to ground state transition, where an 
exciton-photon complex is absorbed into the condensate. 
There is no corresponding excitation in the normal state 
Green's function, because the ground state of the N + 1 
particle system(A~ + 1 excitons) cannot be reached from 
the excited states of the N particle system(iV — 1 excitons 
and 1 polariton) by adding a photon. 

The relationship between the excitation spectrum and 
the phase diagram is slightly different when the transi- 
tion occurs for p ex > 0. For example, in the right panel 



Since the inhomogeneous model has a distribution of 
excitations, we must study the spectral function Aho). 
A(lo) is proportional to the imaginary part of the re- 
tarded Green's function, 



A{uj) = 2ZG R (-uj + p ex ). 



(30) 



It is proportional to the optical absorption coefficient of 
the cavity at frequency to, i.e. the imaginary part of the 
dielectric susceptibility. 

To obtain Aho) we require the retarded Green's func- 
tion G R . In the normal state this is given by the straight- 
forward analytical continuation 



G ho ) = lim Q(iu) n — to — irj). 



(31) 



However, in the condensed state we face the problem, al- 
ready mentioned for the homogeneous case, of separating 
the order parameter response from the quasiparticle re- 
sponse. We circumvent this problem by simply assuming 
that the continuation (|3l|) of the normal state Green's 
function also applies in the condensed state. 

Inverting the Q^ 1 contained in ( f22| ) and using (30) and 
( pl| ) expresses Aho) in terms of integrals over the distri- 
bution of energies of the two-level oscillators. We evalu- 
ate these integrals in the limit 77 — > by setting r\ = in 
the integrands and deforming the contour of integration 
around the poles of the integrand on the real axis. The 
contribution to the integrals from the detour around the 
poles can be performed analytically, leaving a principal 
value integral which we evaluate numerically. 

Figure shows the evolution of our calculated absorp- 
tion spectra, Aho), as we increase the density through 
the transition, for g/3 = 2, a = 1, and A = 0. The cor- 
responding chemical potential is marked as the dashed 
line. For the empty cavity, p ex — —0.5, we recover the 
absorption spectrum calculated by Houdre et al.o. Com- 
parison with Fig. O shows that, for these parameters, the 
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FIG. 7. Spectral functions(optical absorption spectra) A(uj), for A = 0, gj3 = 2, a = 1 and chemical potentials 
(ftex — E g )/g — —5, —1.5, —1.0, —0.76, —0.75, —0.70, increasing from top left to bottom right through the transition at 
(/J-ex — E g )/g — —0.76. The top row of plots are in the normal state, the bottom left hand plot at the transition and the 
remaining plots in the condensed state. The vertical dashed lines mark the chemical potential. 



positions of the polariton peaks are largely unaffected by 
the inhomogeneous broadening. However, since the po- 
laritons are now resonant with a significant density of 
electronic states they become broadened. Increasing the 
chemical potential, but remaining in the normal state, 
we see the thermal occupation factors producing gain 
below the chemical potential and increased absorption 
just above. The collapse of the polariton splitting evi- 
dent in Fig. ^ is hardly noticeable at these low densities. 
As the density is increased still further a pole appears 
in A(u) at the chemical potential; this marks the onset 
of condensation. Above the critical density the coher- 
ent cavity field, oscillating at frequency [i ex , produces a 
gap of magnitude 4g|A| in the spectrum. The peak on 
the high energy side of the gap connects smoothly to the 
upper polariton peak of the normal state, just as in the 
homogeneous case. In the homogeneous case we noted 
the appearance of an excitation below the chemical po- 
tential as we crossed the transition. This is still present 
in the inhomogeneous case, but for the parameters used 
in Fig. it is far too weak to be visible. 

VIII. CONCLUSIONS 

Real microcavities are far more complex than the ide- 
alized model (0). However, like our model, they consist 
of photons coupled to electronic excitations which are 



bosons at low densities, but reveal their fermionic in- 
ternal structure at high densities. We have shown how 
the polariton condensate may be generalized to allow for 
the saturation nonlinearity produced by such fermionic 
structure. The saturation nonlinearity can produce (l)a 
collapse of the splitting between the peaks in the absorp- 
tion spectrum of the normal state with increasing den- 
sity, (2) a shift of the chemical potential of the condensate 
away from the conventional polariton energy, and (3)an 
unusual reentrant phase boundary for condensation. 

Experimental work on cavity polaritons has concen- 
trated on microcavities containing high-quality GaAs 
quantum wells. In these systems, the excitons are 
weakly-bound, and rather delocalized. Thus, while the 
saturation nonlinearity discussed here is present for these 
excitations, it will be accompanied by other nonlineari- 
ties produced by the overlap of the wavefunclims of dif- 
ferent excitons and the ionization of excitonsOEil. These 
effects may well prevent condensation, but are separated 
from the saturation nonlinearity considered here in sys- 
tems with localized, tightly-bound excitons. Note also 
that tightly-bound excitons have a large dipole coupling 
g, and hence the transition temperature will be larger. 

For real examples of localized oscillators, there will be 
some energy E m above which delocalized states exist. 
The picture of a condensate formed from localized oscil- 
lators then only holds when E m — fi ex is large compared 
with /3 _1 and g. By considering Fig. pi we deduce that 
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to completely realize a reentrant phase diagram like that 
shown in Fig. |] requires an energy gap AE separating the 
localized and delocalized excitations; this gap must be 
large compared with g andjS -1 . Such a gap could occur 
in organic semiconductorsQ'Ll. In these systems, excitons 
are strongly bound and therefore small (Frenkel). They 
readily self-trap on local lattice distortions and on im- 
purities in these, often highly disordered, materials. An 
energy gap AE could exist in inorganic quantum wells 
if the excitons move in a potential containing deep, well- 
separated traps, perhaps .asspeiated with interface islands 
in narrow quantum wellsE3~E3. 

Hie disordered quantum wells studied by Hegarty et 
alxH provide an example of a system without a gap sep- 
arating the localized and delocalized excitations. These 
systems show a single inhomogeneously broadened exci- 
ton line, unlike the quantum wells of Refs. |37] [H The 
"mobility edge" E m lies near to the center of the exci- 
ton line. One may be able to form a condensate which 
does not involve delocalized excitations using this type of 
quantum well if the inhomogeneous broadening is large 
compared with g and (3 and the cavity mode is placed 
low down in the exciton line. The transition would then 
occur when the chemical potential is well separated from 
E m . 

The polariton condensate described here is formed 
from a quasi-thermal population of electronic excita- 
tions which are renormalized by the coherent photons 
in the cavity. This renormalization, embodied in the Bo- 
golubov transformation (fl5|), produces a gap of magni- 
tude 4g|A| in the absorption spectrum of the condensate. 
Such renormalizations, and hence ihe gap, are absent in 
conventional semiconductor lasersEj, for which a quasi- 
thermal population of the bare electronic excitations is 
assumed. Thus the presence or absence of a gap allows 
the polariton condensate to be distinguished from a con- 
ventional laser. 

In a conventional laser, the renormalization of the elec- 
tronic excitations by the photons, and hence the gap, is 
absent because the electronic polarization is very heavily 
damped. The destruction of a gap by damping is well- 
known in superconductors, where it is associated with 
magnetic impurities. Such impurities suppress the gap, 
eventually to zero. The destruction of the gap does not 
coincide with the destruction of the order parameter how- 
ever: near to T c there is a regime of gapless superconduc- 
tivity, in which there is an order parameter but no gap in 
the single-particle spectrum. This regime should corre- 
spond to the conventional semiconductor laser, although 
the actual damping mechanisms-, will differ. 

The non-equilibrium analogEa of the crossover illus- 
trated in Fig. |^ and Fig. ^, from a two-peaked polari- 
ton spectrum Id a "Stark triplet", has been observed 
experimentally^. In that experiment, the gapped ab- 
sorption spectrum is observed simultaneously with the 
excitation pulse. Thus there is no thermalization in- 
volved in producing the gapped spectrum. It is the result 
of coherence in the excitation pulse, rather than the spon- 



taneous coherence of condensation. Nonetheless, these 
experiments demonstrate the renormalization of the elec- 
tronic states that is essential in the polariton condensate. 

To reach the quasi-equilibrium regime we have de- 
scribed requires a system where the polariton lifetime 
is long compared with the time required to reach ther- 
mal equilibrium at a fixed number of polaritons. Cur- 
rent semiconductor microcavities have lifetimes for the 
photons, and hence the polaritons, of the order of pi- 
coseconds. Finding an exciton system which thermalizes 
on this timescale may be difficult. However, there seems 
no reason to suppose it is impossible, particularly be- 
yond the linear regime, since nonlincarities can enhance 
relaxationEJ. Furthermore, microcavities are available 
with lifetimes far greater than picoseconds. For example, 
silica microspheres have confined modes with lifetimes of 
microsecondafj. 
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APPENDIX A: THE PHASE MODE 

In this appendix we investigate the zero eigenvalue of 
Q^ 1 that appeared while studying the stability of the 
condensate in the homogeneous case. We first prove that 
the zero is also present in the inhomogeneous model, and 
that it describes phase fluctuations of the condensate. It 
is thus the Goldstone mode reflecting the degeneracy of 
the ground state with respect to the phase of the order 
parameter. We then argue that the zero eigenvalue does 
not contribute to the free energy density in the thermo- 
dynamic limit. Although the physics we discuss in this 
appendix is well understood in general, it is particularly 
transparent in our simple model. 

We note that, at to = 0, K\ is real and positive. 
The eigenvalues of Q^ 1 are then K\ ± \K2\- Since from 
the explicit forms of K\, and the extremal equa- 
tion ( |l6| ) we have \K%\ — Ku as .required in general by 
the Hugenholtz-Pines relatiorcBcj, Q~ x has a zero eigen- 
value. 

To illustrate that the zero eigenvalue is the phase mode 
of the condensate, note that since argi^2 = 2argA = 2cj) 
we can write 



The eigenvector of this matrix with zero eigenvalue is per- 
pendicular in the complex plane to the order parameter 
A. 

Since we are considering a broken symmetry system, 
we should not include states with different phases of the 
order parameter when calculating the partition function. 
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Thus on physical grounds, we should discard the zero 
mode when computing the partition function. 

A formal approach which allows calculations in the 
presence of the zero eigenvalue is to introduce symmetry 
breaking terms which are taken to zero after the thermo- 
dynamic limit. This is the standard method for apply- 
ing statistical mechanics to broken symmetry systemsca. 
The appropriate symmetry breaking terms for a Bose 
condensed system pin the phase of the order parameter. 
They are sources and sinks for the photons, and appear in 
the effective action S c g as -j^ + Jijj) . These terms 

do not contribute directly to (|22|) , but appear as a source 
term in ©. The original zero eigenvalue of Q is now 
K\ — | if 2 1 = —J/(i/>qVN). Since for the equilibrium so- 
lution we must have <p — arg J = 7r, the contribution of 
the original zero eigenvalue to the free energy density is 
proportional to 

lim lim — In ( — ) = 0. 

j^on^ooN \\\\^/nJ 
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